A characterization of horizontal visibility graphs and combinatorics on words 
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An Horizontal Visibility Graph (for short, HVG) is defined in association with an ordered set 
of non-negative reals. HVGs realize a methodology in the analysis of time series, their degree 
distribution being a good discriminator between randomness and chaos [B. Luque, et al, Phys. 
Rev. E 80 (2009), 046103]. We prove that a graph is an HVG if and only if outerplanar and has a 
Hamilton path. Therefore, an HVG is a noncrossing graph, as defined in algebraic combinatorics [P. 
Flajolet and M. Noy, Discrete Math., 204 (1999) 203-229]. Our characterization of HVGs implies 
a linear time recognition algorithm. Treating ordered sets as words, we characterize subfamilies of 
HVGs highlighting various connections with combinatorial statistics and introducing the notion of 
a visible pair. With this technique we determine asymptotically the average number of edges of 
HVGs. 



I. INTRODUCTION 



A graph is an ordered pair G = (V,E), where V is a 
set of elements called vertices and E C V x V — {{i,i}: 
i G V} is set of unordered pairs of vertices called edges. 
Let X — (xi 6 R-° :i = 1,2, ...,n) be an ordered set 
(or, equivalently, a sequence) of non-negative real num- 
bers. The horizontal visibility graph (for short, HVG) 
[13] of X is the graph G = (V, E), with V = X and an 



edge XiXj if 



> Xk for every i < k < j . Clearly 



XiXj 6 E whenever j = i + 1. If the elements of 
X are given with finite precision, we can always take 
X=(ij£N:i = l,2, n). The HVG of X is also de- 
noted by HVG (AT). 

The term HVG also describes some diagrams employed 
in fioorplanning and channel routing of integrated cir- 
cuits [6j ; visibility graphs obtained from polygons consti- 
tute an area of extensive study in computational geome- 
try (see, e.g., Ch. 15 in [l[). Inspired by the geometric 
notion, Lacasa et al. [?J introduced visibilitygraphs and 
shortly afterwards their variations, HVGs [12| ■ Further 
research on such graphs is reported in [1, Q . The princi- 
pal idea of these papers is to translate dynamical prop- 
erties of time series into structural features of graphs. It 
was shown in [12J that the degree distribution of HVGs 
has a role in discriminating between random and chaotic 
series obtained from discrete dynamical systems. In 
the authors use an algorithm for HVGs to characterize 
and distinguish between correlated stochastic, uncorre- 
cted and chaotic processes. 

The purpose of the present paper is to study combina- 
torial properties of HVGs. We characterize the family of 
HVGs by proving that a graph is an HVG if and only if 
is outerplanar and has a Hamilton path. This result al- 
lows us to use, in the study of HVGs, results obtained for 
outerplanar graphs. Recognize properties of a sequence 



X via the properties of its HVG is a wide direction of 
analysis. Treating ordered sets as words, we character- 
ize subfamilies of HVGs highlighting various connections 
with combinatorial statistics. 



II. OUTERPLANARITY 

A drawing of a graph G is a function / that maps 
each vertex v £ V(G) to a point f(v) S M 2 and each 
edge uv 6 E(G) to a curve whose endpoints are f(u) and 
f(v). A graph is planar if it has a drawing so that any 
pair of edges can only intersect at their endpoints. A face 
is a bounded region of a planar graph. The infinite face 
is the outer one. A planar graph is outerplanar if it has 
a drawing such that each vertex is incident to the infinite 
face. A linear time recognition algorithm for outerpla- 
nar graphs has been given by Mitchell [l4| . Outerplanar 
graphs are used in the design of integrated circuits when 
the terminals are on the periphery of the chip [101 ] . In this 
section, we give a characterization of HVGs with respect 
to outerplanarity. The characterization is a consequence 
of the following two results: 

Lemma 1 Every HVG is outerplanar and has a Hamil- 
ton path. 

Proof. Let G = HVG(AT), where X = ( Xl ,...,x n ). 
By definition, G has a Hamilton path P = Xix% . . .x n . 
Observe that there is no a pair XiXj, x s Xt of edges such 
that i < s < j < t. Indeed, the existence of such a pair 
implies that Xj > x s (due to XiXj S G) and Xj < x s (due 
to x s Xt G G), a contradiction. Thus, we can obtain an 
outerplanar embedding of G as follows: set as an interval 
of a horizontal straight line with points being vertices of 
P and depict each edge of G outside P as an arc above 
the interval for P such that if x%Xj, x s xt € E{G) \ E(P) 
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and s < i < j < t then the arc corresponding to XiXj is 
situated below the arc corresponding to x s x t . Hence, we 
have proved that G is outerplanar. ■ 

Lemma 2 If a graph H is outerplanar and has a Hamil- 
ton path, then H is an HVG. 

Proof. Let H have a Hamilton path V\v^ ■ ■ ■ V n . As H is 
outerplanar, all vertices of the Hamilton path belong to 
the infinite face in an outerplanar embedding of H into 
plane. Observe that H has no crossing edges, i.e., pairs 
of edges ViVj, v s v t such that i < s < j < t. We say that 
an edge ViVj, i < j, of H is of nesticity if j = i + 1. 
We say that an edge ViVj, i < j, of H is of nesticity 
k if k is the minimum nonnegative integer such that if 
v s v t € E(H) and i < s < t < j then v s v t is of nesticity 
p, where p < k. Nesticity of each edge is well-defined as 
there are no crossing edges. The nesticity of a vertex Vi 
is the maximum nesticity of an edge incident to m. 

Now we will construct an ordered set X — (xi, . . . , x n ) 
such that H = HVG(V) using the following algorithm 
which consists of two stages. In the first stage, we deter- 
mine the nesticity of each edge of H . To initialize, set the 
nesticity of each edge XiXi+i to 0. Now for each q from 
2 to n — 1 consider all edges of H of the form XiXi+ q and 
set the nesticity of XiXi+ q to p + 1, where p is the max- 
imum nesticity of an edge XjXk with i<j<k<i+q 
and XjXk ^ XiXi +q . In the second stage, compute the 
nesticity of each vertex m of H by considering all edges 
incident to V{ and set Xi to the nesticity of Vi. 

It is easy to see, by induction on the nesticity of edges, 
that indeed H = HVG(V) (we use the fact that H has 
no crossing edges). ■ 

Theorem 3 A graph is an HVG if and only if it is out- 
erplanar and has a Hamilton path. 

Since recognizing an outerplanar graph and determin- 
ing if it has a Hamilton path are tasks that can be done 
in linear time (see Mitchell 141] and Lingas [llj, respec- 
tively), we have following result: 

Corollary 4 We can determine in linear time (with re- 
spect to the number of vertices) if a graph is an HVG. 

We conclude this section with a further characteriza- 
tion. A noncrossing graph [5[ with n vertices is a graph 
drawn on n points numbered in counter-clockwise order 
on a circle such that the edges lie entirely within the cir- 
cle and do not cross each other. It is useful to point out 
that these objects are HVGs because of Theorem [3] 

Corollary 5 An HVG is a noncrossing graph. 

III. UNIMODAL HVGS 

We characterize here the HVGs with minimum num- 
ber of edges. The degree of a vertex i is d(i) := 



\{j '■ {^j} £ The degree sequence is the unordered 

multiset of the degrees. HVGs are not characterized by 
their degree sequence. In other words, there are non- 
isomorphic HVGs with the same degree sequence. For 
example, given X — (xi : 1, 2, 5), let x\ — X3, X2 < x\ 
and x 2 — X4 = x 5 . The degree sequence of HVG(A) 
is {3,2,3,2,2}. The same degree sequence is associated 
to the graph HVG(A'), where x[ = x' 4 , x' 2 < x\ and 
x 2 = x 3 = x 5- Let (5(G) and A(G) be the minimum 
degree and maximum degree of a graph G, respectively. 
Theorem [3] implies the following: 

Proposition 6 If G is an HVG then 5(G) = 1 or 2 and 
A(G) < n— 1. If A(G) = n — 1, then there is only one 
vertex of maximum degree. 

A real-valued function / is unimodal if there exists a 
value m such that / (x) is monotonically increasing for 
x < m and monotonically decreasing for x > m. Thus, 
/ (ra) is the maximum of / and its only local maximum. 
If X is unimodal then HVG(X) is said to be unimodal. 

Proposition 7 An HVG is unimodal if and only if it is 
a path. 

We associate a word of length n — 1, called a differ- 
ence, to a set with n elements. The letters of the word 
are from the alphabet {0, +, — }. Each letter corresponds 
to a pair of adjacent elements. If Xi = ajj+i the letter cor- 
responding to the pair (2^,2^+1) is 0. When Xi < Xi+i 
(resp. Xi < Xi+i) then the letter for the pair is -I- (resp. 
— ). For example, X = (1, 4, 8, 8, 2, 4) gives the difference 
+ + — h Let us observe that the difference D of a 
unimodal HVG graph does contain the pattern (—,+), 
i.e., it does not contain a subword — P+, where P is an 
arbitrary subword of D. 

Proposition 8 The number of HVGs isomorphic to the 
n-path and associated to different ordered sets of cardi- 
nality n ( without taking into account the actual values of 
the single elements) is exactly a(n) = 2 n ~ 1 (n + 2). 

Proof. We prove that the number of differences of 
length n without the pattern (— ,+) is exactly a(n). 
Clearly, oq — 1 and a\ = 3. Since each such a word 
x = X1X2 ■ • ■ x n of length n can be written as x = +x' , 
x = Ox', or x — —x" then a n = 2a„_i + b n , where b n 
denotes the number words in B n of length n on the al- 
phabet {0, +, — } that do not contain the pattern (— , +) 
and its leftmost letter is — . Since any word in B n has 
leftmost letter — and there is no pattern (— , +) in the 
word, any letter which is not the leftmost has two possi- 
bilities: or — . This implies that b n — 2 n ~ 1 and, hence, 
a n satisfies the recurrence relation a n — 1a n —\ + 2™" 1 
with initial conditions ao = 1 and a\ = 3. Solving this 
recurrence relation gives the formula for a n . ■ 

The difference does not uniquely specify the HVG. In 
fact, there are nonisomorphic HVGs with the same dif- 
ference. For example, the sets (5,4,3,5) and (5,4,3,4) 

have the same difference h, but the HVGs associated 

to these sets have five and four edges, respectively. 
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IV. MAXIMAL HVGS 

The triangulation of a (convex) polygon is a planar 
graph obtained by partitioning the polygon into disjoint 
triangles such that the vertices of the triangles are chosen 
from the vertices of the polygon. An outerplanar graph is 
maximal outerplanar if it is not possible to add an edge 
such that the resulting graph is still outerplanar. A max- 
imal outerplanar graph can be viewed as a triangulation 
of a polygon. By Theorem 

Corollary 9 The maximum number of edges in an HVG 
on n vertices is 2n — 3 and the graph is maximal outer- 
planar. 

In computational geometry, the visibility graph (TEl ] of 
a polygon with n angles is obtained by constructing a 
graph on n vertices, each vertex of the graph represent- 
ing an angle of the polygon, and each edge of the graph 
joining only those pairs of vertices that represent visible 
pairs of angles in the polygon. A polygon in the plane 
is called monotone with respect to a straight line L, if 
every of the lines orthogonal to L intersects the polygon 
at most twice. By a result of ElGindy [4] and Theorem[3j 
we have an observation relating HVGs to a special classes 
of visibility graphs: 

Corollary 10 An HVG on n vertices and 2n — 3 edges 
is the visibility graph of a monotone polygon. 

A characterization of HVGs with maximal number of 
edges is given with Corollary Q1JJ In fact, the charac- 
terization is easy once established a connection between 
HVGs and combinatorics on words. 



Theorem 12 Let X — {xi, X2, ■ ■ ■ , x n } be an ordered 
set of n elements and let k = max A. The HVG of 
X , G =HVG[X) , can be represented uniquely as a word 
x = x\X2 ■ ■ ■ x n with a set E of visible pairs, where each 
edge {i,j} in the graph G corresponds to visible pair 

(Xi ,Xj). 

Let A = (xi, Xj ) and B = (xi> , x.y ) be two visible pairs 
in a word x — x\X2 • • • x n £ [k] n with i < i' . Clearly, 
i < j and i' < j'. From the above definition, we obtain 
that either i < i' < j' < j or i < j < i' < j' . In such 
cases, we say that the pair A covers the pair B and that 
A and B are disjoint pairs, respectively. 

Fact 13 If x is the word realizing an HVG then in any 
two disjoint pairs A, B £ x either A covers B or B covers 
A. 

How many words are there in [k] n with a fixed number 
of visible pairs? 

Definition 14 We denote the generating function for 
the number of words x £ [k] n according to the number 
of visible pairs in x by F}.{x, q), that is, 

F k {x,q) = Y, xn E 

Similarly, we denote the generating function for the num- 
ber of words x £ [k—l] n according to the number of visible 
pairs in xk (respectively, kx) by L k (x,q) (respectively, 
Rk(x,q)). Also, we denote the generating function for 
the number of words x £ [k— 1]™ according to the number 
of visible pairs in kxk, which are not equal to (k,k), by 
M k (x,q). 



Note that each word x in [k] n can be decomposed either 



as 



V. WORDS 

We denote by [k) n the set of all words of length 
n over an alphabet [k] = {1,2, ... ,k}. Each word 
x = x\X2 ■ ■ ■ x n defines an ordered set X by X = 
{xi, X2, ■ ■ ■ , x n } and conversely. In order to describe the 
edges of an HVG with respect to words, we need the 
following definition: 

Definition 11 Let x — x\X2 ■ • ■ x n be any word in [k] n . 
We say that the pair (xi , Xj ) with i + 1 < j is visible if 
x i+ i,x i+2 , ■ ■ ■ ,Xj-i < min{xi,Xj}. Clearly, (xi,x l+ i) is 
a visible pair. We denote the number of visible pairs in 
x by vis{x). 

For instance, if x = 21232143112112 is a word in [4] 14 
then (2:1,2:3), (2:4,2:7), (x5,X7), (xs,Xii) and (2:11,2:14) 
are the visible pairs of x, in addition to the 13 edges of 
P14. Thus vis(x) = 18. 

By the definition of HVG and visibility of pairs in 
words we can state the following result. 



• x £ [k — 1]", that is, x does not contain the letter 

k, 

• x = x^kx^ ■ ■ ■ kx^ m+1 \ where x^> is a word in 
[k — with £1 + • • • + im+i + m = n. 

Thus, rewriting these rules in terms of generating func- 
tions we obtain the following lemma. 



Lemma 15 For all k > 1, 

Fk(x, q) = F k -i{x, q) + 
and for all k > 2, 



M k (x,q) = M k -i(x,q) + 



xRkjx, q)L k {x, q) 
l-xqM k (x,q) 

xq 2 M k _ 1 (x,q) 



1 - xqM k -i(x,q) 



n / \ R k -i(x,q) 

Rk(x,q) = — r 

1 - xqM k ^i(x,q) 

L k (x,q) = R k (x,q). 
Also Fq(x, q) = Mi(x, q) = R\ {x, q) = L\{x, q) = 1. 
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The above lemma together with induction give the fol- 
lowing result. 

Theorem 16 Let k > 1. The generating function 
F k (x, q) is given by 



F k (x 



l-xqMj(x,q) 



-xqMi(x,q)) 2 
where M k (x,q) satisfies the recurrence relation 

xq 2 M k _ 1 (x,q) 



M k (x,q) = M fc _i(ar,g) + 



1 - xqMk-i(x,q)' 
with initial condition Mi(x,q) = 1. 
For instance, the theorem gives 

F 1 (x,q) = l 

and 

F 2 (x,q) = 1 



1 — xq 
x((l - xq) 2 + 1 -x 2 q 3 ) 



(1 — xq)((l — xq) 2 — x 2 q 3 ) ' 

Note that it appears to be hard to derive an explicit for- 
mula for the generating functions M k (x,q) and F k (x,q). 
But we can use the result for studying the total number 
of visible pairs in all words in [k] n . 

Theorem 17 The generating function for the total num- 
ber of visible pairs in all words in [k] n (k > 1) is 

Fi(x) = j2* n E vi < x ) 

n>0 xe[k] n 



A=l (X-{i-X)x){l—ix) 

j= - (l-{j-l)x)(l-jx) 
k 1 (n l)x I x 2 -( 2l ~ 1 ) x 

9\T^_ _ ^^=1 l-(t-l)a 

x 4- a 



£ (l-(j- 1)1)2(1-^)2 



Proof. We define 



and 



M' k (x) = —M k (x,q) |, =1 



Theorem 1 1 61 with induction on fc implies that F k (x,l) — 
— ^ and M k (x,l) = 1 _ (k 1 _ 1)x - By differentiating re- 
spect to q and by using the expressions F k (x,l) and 
M k (x, 1), Theorem [TBI gives 



^(x) = -x 2 V — -f 

fen J *=i(i 



Mi + Mj 



2l 2f \-xMj J^ Mj + Mj 



J= inLi(i-^) 2 tt 



where Mj — 1 _^ 1 _ 1 ^ x and Mj is given by the recurrence 
relation 

_ (l-(fc-2»2 
Aife ~ (l-(k-l)x)* Mk - 1 



+ 



(i-(k-i) x y 

x(2 - (3k - 3)x) 
(1 - (k-2)x)(l - (k - l)x) ' 



with the initial condition M[ = 0. Hence, by induction 
on k we obtain that 



ML 



^ 2 - (2j - l)x 



(1 - (k - l)x) 2 1)= 



Plugging this into the equation of F k we have the result. 
For instance, ■ 

• F{(x) = 

- Z?'fV\ — (4-3x)x 2 

• *2\ x ) ~ (l-x){\-2x)*' 



H{x) 



{10x 2 -22x+9)x* 
(1-x)(1-2k)(1-3k) 2 



Theorem [18] shows that the function F k (x) has a pole 
of degree two at jr. Direct calculations show that 

C= \im{l-kx) 2 F' k {x) 



lim 



x 2 (i-ik-^x + xYHzl 2 -^^ 



-(i-l)x 



(l-(k-l)x) 2 



1 ^ 2k + 1 - 2i 
k ^ k+l-i 

i—l 

tf(fc + l) +7 



2- 



where \l/(a;) is the digamma function and 7 is Eulcr's 
constant. Thus, asymptotically, the total number of vis- 
ible pairs in all words in [k] n (k > 1) is given by Gnk n . 
This means that the average number of visible pairs in 
all words in [k] n is given by Cn. 

Corollary 18 The average number of edges in an HVG 
is 

(2 -(#(*+ 1) + 7)/*) n, 

when X is an ordered subset of n elements from [k] and 
n — ^ 00. 

The corollary shows that the number of maximal edges 
is in fact 2n — 3, where X contains n elements (cfr. Sec- 
tion III). This can be obtained easily from our represen- 
tation as words: 

Corollary 19 Let G be an HVG with 2n — 3 edges. Then 
X = {..., 8,6, 4, 2, 1,3,5, 6,7,9,...} 

or 

X = {... ,9,7,5,3, 1,2,4,6,8,...}. 



5 



Proof. Without loss of generality, we can assume that 
X has n different numbers. Assume tt = ttitt2 ■ ■ • tt„ is 
a permutation on [n] with a maximal number of visible 
pairs. The pair (tti, tt u ) is visible, hence cither tt\ + 1 — 
TT n = n or 7r„ + 1 = it i = n. Delete n from tt and 
denote the resulting permutation by tt' By induction on 
n, we obtain that tt has the form tt = ■ ■ ■ 864213579 • • ■ 
or tt = ■ ■ ■ 975312468 • ■ ■ , as claimed. ■ 



VI. PERMUTATIONS 

We denote by S n the set of all permutations on [n\. 
Each permutation tt = tt\tt2 ■ ■ ■ 7r„ defines an ordered set 
X by X — {7Ti, 7T2, . . . , 7r„} and conversely. By definition 
im we can describe the edges of an HVG with respect to 
permutations (a permutation is in fact a word without 
repetitions). Define S n (q) to be the generating function 
for the number of permutations tt on [n] according to the 
number of visible pairs in tt by S n (q), that is, 



S n (q) 



For example, in S2 there are two permutations 12 and 
21. Thus S2(q) = 2q. In S3 there are 6 permutations 
123, 132, 213, 231, 312 and 321. Thus S 3 (q) = 4q 2 + 2q 3 . 

Now, let us find an explicit formula for S n (q). Define 
S n ,j to be the set of permutations TT — ^1^2 ' * ' ^"n In Sn 
such that tt j = 1. Let tt 6 S n and define tt' to be the 
permutation obtained from tt by deleting the letter 1 and 
by decreasing each letter by 1. Let us write an equation 
for S n (q). From the definitions, we have that, if j = l,n 
then the set permutations S n ,j is counted by qS n -i{q), 
and if 2 < j < n — 1 then the set of permutations S n j is 
counted by q 2 S n -\{q). Thus, for all n > 2, 

S n (q) = 2qS n -i(q) + (n - 2)g 2 5„_i(g), 

which implies that 

Sn (q) = (2qr- 1 f[(l + J —l, 

Using the fact that the unsigned Stirling numbers s(n, j) 
of the first kind satisfy the relation 



(1 + a:) • ■ • (1 + (n - l)x) 



y%(n,n- j)x j , 



we obtain that 



n-2 



S n (q) ^ (2q) n - 1 U (l+J- 
3=1 



■ 1 



n—1 j 

(2g)«- 1 ^ S (n-l,r 1 -l-i)| J 
3=0 



n-l 

E 

3=0 



2 n - 1 - j s(n - l,n-l-j)q 



n—l+j 



Hence, we can state the following result: 

Theorem 20 Let n > 2. The number of permutations 
tt with exactly n—l+j , < j < n — 1, visible pairs is 
given by 

2«- 1 -'s(n-l,n-l-j), 

where s(n — 1, n — 1 — j) is the unsigned Stirling number 
of the first kind. 

Corollary 21 The average number of edges in an HVG 



1 

when X is an ordered subset of n different elements. 

Proof. By Theorem[20]we have that the average number 
of edges in an HVG is 

n-l 

Pn = -. Y> - 1 + j)2 n - 1 - J s(n - l,n - 1 - j), 
3=0 

when X is an ordered subset of n different elements. The 
expression p n can be written as 

n — 1 

Pn = ~. V(2n-2-j)2^ S (n-l,j). 



3=0 



By the fact that 



x(a; + 1) • • • {x + n — 1) = s{n,j)x^, 

3=0 



we obtain 



ft lb — -L 

y~]js(n, j)x J = x(x + 1) • • • (x + n - 1) V — — , 
^-^ ^ — ' x + 1 

3=0 j=0 J 



which implies that 



J2 s ( n - l dW = n! 

3=0 



and 



n-l 



J2j^s(n-l,j) = n\J2-. 



3=0 



3=2 



Hence, 



™ 1 n 1 

Pn = 2(n - 1) - ^ 7 = 2n - E 7' 

which completes the proof. 



3=2 3 3=1 3 
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VII. CONCLUSIONS 

We have characterized the family HVGs in terms of 
their combinatorial properties. The characterization is 
useful because it gives an efficient recognition algorithm. 
Moreover, it places the study of HVGs in a specific math- 
ematical context, related to a well-known class of graphs. 
As it is originally observed in the literature, the poten- 
tial importance of HVGs stems in their use for describing 
properties of dynamical objects like time series. There- 
fore, the main goal would be to determine the dynami- 
cal and structural properties of an ordered set that are 
readable through the analysis of its HVG. In this per- 
spective, we have shown that combinatorics on words are 
a useful tool. The connection suggests a number of nat- 
ural open problems. For example, it may be valuable to 
study HVGs and statistics on words and forbidden subse- 
quences. The mathematical scope turns out to be wider 



than the original dynamical systems framework. 

It would be interesting to characterize visibility graphs 
0, Hj], which can be defined as follows. For an ordered 
set X — (xi G R-° : i = 1,2, its visibility graph 

VG(X) has the vertex set X and Xj, Xj (i < j) are adja- 
cent if for each k with i < k < j we have 

Xfc <C Xj -p \Xj Xj)— 

3-t 

It is easy to see that for each X, HVG(AT) is a subgraph 
of HV(X). Not each visibility graph is outerplanar as 
HV(Y), where Y = (4,2,1,4) is K t , which is not outer- 
planar. 

Acknowledgments. We would like to thank Lucas La- 
casa for useful conversation about HVGs. 
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